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ABSTRACT 

We have investigated three SNRs in the LMC using muhi-wavelength data. These SNRs 
are generaUy fainter than the known sample (see §4) and may represent a previously missed 
population. One of our SNRs is the second LMC remnant analyzed which is larger than any 
Galactic remnant for which a definite size has been established. The analysis of such a large 
remnant contributes to the understanding of the population of highly evolved SNRs. 

We have obtained X-ray images and spectra of three of these recently identified SNRs using 
the XMM-Newton observatory. These data, in conjunction with pre-existing optical emission-line 
images and spectra, were used to determine the physical conditions of the optical- and X-ray- 
emitting gas in the SNRs. We have compared the morphologies of the SNRs in the different 
wavebands. The physical properties of the warm ionized shell were determined from the Ha 
surface brightness and the SNR expansion velocity. The X-ray spectra were fit with a thermal 
plasma model and the physical conditions of the hot gas were derived from the model fits. Finally, 
we have compared our observations with simulations of SNR evolution. 

Subject headings: supernova remnants — X-rays: diffuse background 

1. Introduction for developing a complete understanding of the 

complex structure and evolution of the ISM in 
Supernovae (SNe) and their remnants are the gas-rich galaxies not unduly affected by external 
driving force behind the dynamic interstellar ^^^^^^ ^^^^ ^s galaxy-galaxy collisions.. However, 
medium (ISM) of our Galaxy and others. The order to form an accurate picture of how SNRs 
study of the energetics and evolution of super- influence the ISM, it is necessary to have a com- 
nova remnants (SNRs) thus forms the foundation p^g^g ^nd unbiased sample of SNRs in a galaxy. 
This is very difficult within the Milky Way as dust 

1 Present email address: klimek@physics.utexas.edu absorption and COnfusion along the line-of-sight 
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hinders our ability to detect and determine the 
physical conditions of Galactic SNRs. Although 
it is possible to obtain a more complete sample of 
SNRs in more distant galaxies, their distance hin- 
ders the investigation of their physical conditions 
in detail. 

The Large Magellanic Cloud (LMC) offers the 
ideal laboratory for the study of a large sample 
of SNRs. The LMC is well out of the Galactic 
plane and thus offers low foreground absorption 
{Ay ^ 0.3). Meanwhile, its distance is re l atively 
well k nown (~50 kpc: iDi Benedettol (|2008l ). iFeast 
and distances within the LMC are small 
compared with its distance to us so all objects 
within the LMC can be taken to be at approxi- 
mately the same distance. Its proximity also al- 
lows us to resolve a great amount of spatial detail 
(r'« 0.25 pc). 

The SNRs presented in this paper were discov- 
ered as regions of shock-ionized gas in the Mag- 
ellanic Cloud Emission Line Survey ( MCELS) us- 
ing t he reliable [Sii]/Ha diagnostic ()Fesen et al. 
19851). with SNRs tvpicallv having [Sii]/Ha > 
0.4 ( Levenson et al.l fl995i) . Analysis of archival 



ROSAT data and high resolution optical spec- 



tra confirmed these shells as SNRs (ISmith et al 



(|2011. in preparation! ) . ISmith et al.l (|2003[) V How- 
ever, the ROSAT observations are too shallow to 
provide much detailed information on the hot gas 
within these remnants. Two of these remnants 
have been detected, but not classified as SNRs in 



the ROSAT survey of the LMC (iHaberl fc Pietsch 



1999 1 ) . a nd the Park e s radi o survey (jFilipovic et al 



19981). iBlair et al 



(|2006[ ) seached for FUV em- 
mission from these remnants but failed to detect 
them. 

A thorough investigation of the hot gas is cru- 
cial to understanding the properties of the SNR 
and its interaction with the host galaxy. SNRs 
contribute most of the hot gas component of the 
ISM, and much of the energy originally released in 
the supernova event is carried as thermal energy 
in the hot gas. Moreover, the interaction between 
the hot interior and the cool shell of shocked ma- 
terial dictates the evolution of the SNR. We have 
therefore obtained new X-ray data from the XMM- 
Newton observatory and analyzed these data along 
with the optical data from MCELS to form a com- 
plete picture of these SNRs. 

Our observations and analysis techniques are 



presented in §2. In §3 we give the results of these 
techniques applied to the individual objects, fol- 
lowed by discussion of those results in §4. A brief 
conclusion is given in §5. 

2. Observations 
2.1. X-ray 

XMM-Newton is an orbiting X-ray observa- 
tory operated by the European Space Agency. In 
this study, we use data from the European Pho- 
ton Imaging Camera (EPIC) MOS CCDs. These 
are two identical devices each mounted behind 58 
nested Wolter type 1 mirrors. The EPIC cameras 
offer the capability to perform extremely sensitive 
imaging observations over the telescope's field of 
view of 30 arcmin and in the energy range from 
0.15 to 15 keV with moderate spectral {E/AE ^ 
20 — 50) and angular resolution (point-spread func- 
tion: 6 arcs ec FWHM). For fur ther technical dis- 
cussion, see [Turner et al. I (1200 ll) . 

This study presents new X-ray data from 
XMM-Newton on three LMC SNR candidates: 
SNR0449-6921, SNR0506-6542, and SNR0537- 
6628. Proposed pointings were chosen to follow 
up on regions of high [Sii]/IIa from the MCELS 
data. Three observations are presented in this 
paper. See Tabled] for details on the observations. 

These data were processed with the Extended 
Source Analysis Software (ESAS) package. ESAS 
is a package of perl scripts which call commands 
from ESA's Science Analysis Software and stan- 
dalone FORTRAN 77 progr ams following ana l- 
ysis techniques given by .Snowden et al.l (|2008l ). 



ESAS automatically filters times of high back- 
ground contamination and poor event grades. 

The low X-ray surface brightness of these SNRs, 
in combination with limits of X-ray coverage and 
resolution from past surveys, explains their non- 
detection until the present. We will show that the 
average X-ray surface brightness of these remnants 
falls below the lowest levels detected from LMC 
SNRs included in earlier catalogs. Given the rela- 
tively small number of source counts we are able to 
collect from the remnants, we need to have a good 
understanding of the backgrounds so that we can 
accurately extract source properties from limited 
data. 

After the data had been filtered, we applied the 
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ESAS cheese task, which uses SAS tools to iden- 
tify point sources, and creates a mask file. We 
specifically ran the source detection over a band 
from 0.8 to 8.0 keV in order to identify trou- 
blesome nonthermal background sources without 
falsely masking softer lumpy features within the 
SNR emission. 

Using these masks, we then extracted spectra 
from background regions near the remnants. We 
binned these background spectra to 50 counts per 
channel and fit them with NASA's Xspec package 
using a model composed of two gaussian compo- 
nents to account for the strong Al and Si Ka fiu- 
orescent lines near 1.5 and 1.75 keV, respectively 
an APEC component for the Local Hot Bub- 
ble/Solar Wind Charge Exchange, an absorbed 
APEC component for the Milky Way's halo, and 
an absorbed power law component for the non- 
thermal extragalactic background. The photon 
index of the extra galactic pow er law w as fixed at 
1.46 (jSnowden et al. 2008 : C hen et all 17997.) . Fi- 
nally, we also included a power law component 
which was not folded with the instrumental re- 
sponse to model residual soft proton contamina- 
tion not removed in the light curve filtering. The 
specifics of the background modeling are discussed 
below for each observation, but aside from the soft 
proton contamination which varies with time, the 
background parameters found by this method are 
consistent between all the observations. In addi- 
tion, our temperature for the local hot bubble in 
SNR0449-6921 an d SNR0537-6628 ar e consistent 
with the resuhs of iHenley k SheltonI HoOS). An 
example of our background model is seen in Figure 

m 

Next we extracted spectra of the SNR emission 
and fit each with the corresponding background 
model, plus an APEC thermal plasma component, 
with abundances set to the canonical LMC value 
of 0.3 solar. We froze the temperatures, photon 
indices, column density and gaussian parameters 
as determined from the background fit, but al- 
lowed the normalizations to vary to account for 
the change in area between the background and 
source spectrum extraction regions. The source 
spectra generally have few energy bins with signif- 
icant numbers of counts, but the constraints pro- 
vided by the background fits allowed us to acheive 
statistically significant fits for the source parame- 
ters. Fits using the older Raymond-Smith model 



or the MEKAL model in place of the APEC model 
gave nearly identical results. After this initial fit- 
ting, we added the additional constraint that the 
ratio of the normalizations of the absorbed ther- 
mal and nonthermal components should remain 
constant over the field-of-view. This constraint 
did not change the best-fit values that we found, 
but did improve the uncertainties in the fitted val- 
ues. An example of the residuals from our back- 
ground model fitted with the APEC model is seen 
in Figured 

We compared our fitted values of hydro- 
gen column density with that found in H i sur- 
veys. We used the Leiden/Arge ntine/Bonn sur - 



vey of Galactic hydrogen (|Kalberla et al.l l2005l ) 



and the ATCA-|-Parkes survey of LMC hy- 
drogen (Kim et al. '2003V It has been shown 
(jArabad jis fc Bre gman 1999) that above Nn > 
5 X 10^*^ cm~^, the total X-ray absorption col- 
umn increases by about a factor of 2 over the 
H i-derived value due to molecular hydrogen. We 
must also remember that only part of the LMC 
column is actually in front of our objects. Thus 
we should expect our fitted values to be close to 
2ArjjGai ^ ^jjLMC^ where N^^ is the value quoted 

from the respective survey. In every case, we found 
that our fitted values of A^h are consistent with 
the survey values. 

Finally we compared the observed X-ray spec- 
trum with simulations. We simulated each SNR 
with the aid of a computer code that followed 
the hydrodynamics of the plasma on a spherically- 
symmetric Lagrangian mesh, simulated the shock 
front, shock heating, thermal conduction of heat 
from the remnant's hot center outwards, radiative 
cooling by non-collisional ionizational equilibrium 
(non-CIE) plasma, and, when needed, nonther- 
mal pressure due to magnetic fields and cosmic 
rays. The products of our simulations, the tem- 
perature, density, and ionization level as a func- 
tion of location within the SNR and SNR age, 
were used to calculate the spectra in a non-CIE 
manner . The a lgorithms are described in detail in 
Shelton (19981). The spectra were then converted 



to Xspec models with age as a free parameter. To 
accomplish the spectral fitting in Xspec, we used 
the fit to the source region as described above, 
and assumed the APEC source component accu- 
rately models the SNR plasma emission. We then 
deleted the APEC component and replaced it with 
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the simulated spectrum. We stepped through the 
spectra at various ages and searched for a mini- 
mum in the Chi-squared. The simulations were 
run assuming solar abundances. As a first ap- 
proximation, we varied the normalization of the 
resulting spectral model to account for the differ- 
ent abundances of the LMC. 

We also used ESAS to produce combined and 
adaptively smoothed images in three bands (0.35- 
0.85, 0.85-2.0, and 2.0-8.0 keV) with a kernel of 
15 counts. The ESAS task proton uses the soft 
proton contamination model parameters found in 
the background fitting to subtract this component 
from the images. 

Given the thermal plasma model normalization 
A we can calculate the electron density in the hot 
gas by Ue = 3.89 x lO'^D^/A/Vhot where D is the 
distance to the SNR, A is the normalization of 
the plasma model component, and Vjiot is the vol- 
ume of the hot gas in cgs units. The mass of the 
hot gas is Mhot = l-17nempVhot, assuming fully 
ionized helium and where rrip is the proton mass. 
Given the temperature from our spectral fits, the 
total thermal energy is then E — {3/2)NkT = 
4.60 X lO^^TieVhotT, where T is the temperature 
of the plasma model in keV. We can also find the 
pressure as Phot = nkT = 3.05 x lO'^ngT. These 
calculated values are listed in Table [21 The er- 
rors are the one sigma uncertainties estimated by 
Xspec in the fit process propagated through the 
above equations. 

2.2. Optical 

Optical images c ame from the MCELS project 



sky. The spatial coverage along the slit is roughly 
3', limited by the optics of the camera. Details 
of the reduction of these data can be found in 



Smith et al.l (|2011. in preparatioij) 



( Smith et al. 20001) where these remnants were 
originally identified. The observations consist of 
red and green continuum images and Ha, [Sii], 
and [O III] emission line images from the Cur- 
tis Schmidt Telescope at the Cerro Tololo Inter- 
American Observatory (CTIO). The final images 
are flux calibrated and rescaled to have 2" pixels. 
The continuum exposures were used to subtract 
the continuum contribution to the emission line 
images. 

In addition, we have long-slit, high-dispersion 
echelle spectra of the SNR candidates from the 
CTIO 4m telescope. The data have a 24/Ltm pixel 
size that corresponds to 0.082 A (3.65 km s~^) 
along the dispersion axis and ^ 0'.'26 along the 



The echelle profile of the Ha line is broadened 
along the dispersion axis both by the inherent in- 
strumental profile (~ 14 km s^^) and by thermal 
Doppler broadening ('^ 18 km s~^ at X-ray tem- 
peratures). To assign an expansion velocity to 
each remnant, we extract profiles along the dis- 
persion axis where the emission shows the greatest 
dispersion and measure the peaks of the profile. 

Using the fiux-calibrated MCELS images, we 
measured an Ha surface brightness and the aver- 
age radius R and thickness Ai? of the shell. For 
a uniform spherical shell, the greatest line of sight 
through the shell is £ = 2^ - {R - I\Rf . The 
measured surface brightness then implies an emis- 
sion measure EM = nIC = 5 x 10^'' x SB where 
SB is the surface brightness in cgs units and arc- 
seconds. The total mass of the warm ionized shell 
is M — \.27nf.mpVshc\\i where T4hcii is the volume 
of the shell and assuming singly ionized helium. 
Expansion velocities Woxp have been determined 
from echelle spectra of the Ha emission-line, and 
thus we may determine kinetic energies for the 
warm shells hy K = Mv'^^p/2. If we assume that 
T = 10'' K, we can also calculate the pressure in 
the sheU as Psheii = 2nekT = 2.76 x lO'^^^^^ 

In calculating the volume of the shell, we as- 
sumed an ellipsoidal geometry, measuring two axes 
from the projected face of the shell and taking the 
third line-of-sight axis to be the average of these 
two. We then took the uncertainty in the radius 
along this third axis to be the deviation of the 
first two axes from their mean. We also consid- 
ered the instrumental and Doppler broadening of 
our echelle spectra as an additional source of er- 
ror. Convolving the two profiles, we found an un- 
certainty of 11 km s^^. Calculated quantities are 
listed in Tabled 

The age of the SNR may be estimated us- 
ing the analytic expressions of the Sedov- Taylor 
solution for blast wave expansion, in which the 
shell radius is given as a function of time by 
r{t) = 1.17(£:t2/p)i/5^ where E is the explo- 
sion energy and p is the ambient density. Tak- 
ing the time derivative and applying the Rankine- 
Hugoniot conditions for a strong shock to relate 
the blast wave velocity to the post-shock gas ve- 



4 



locity u, we find u{t) = 0.351{E/pt^y/^. Solving 
ttiis system yields t = 0.3(r/u). We then obtain 
the ages by plugging in the radius observed in the 
optical images and the velocities from the echelle 
spectra. These procedures apply best to remnants 
that are still in the Sedov- Taylor phase of evolu- 
tion, since post-Sedov-phase remnants continue to 
expand at a slower rate. The Sedov- Taylor phase 
is expected to end when the post-shock gas ve- 
locity drops to ~ 190 km s~^. Inhomogeneities in 
the ambient medium also affect the expansion rate. 
In addition, the observed velocities and radii are 
found for individual parts of the SNRs and, due to 
asymmetries in the SNRs may not adequately rep- 
resent the whole objects. All of these effects com- 
bine to make the ages estimated from t = 0.3(r/u) 
approximate. 

The simulations described previously also pro- 
vide an independent numerical estimate of the age. 
Alternatively, given the present values of r and u, 
one may solve for the ratio E/p, or equivalently 
E/n where n is the ambient number density, and 
compare this with the value estimated using the 
X-ray spectral modeling and Ha surface bright- 
ness techniques as described above. 

3. Results 

3.1. SNR0449 6921 

SNR0449-6921 shows the simplest structure of 
the three remnants presented here. Its well defined 
[S ii] shell surrounds a smooth elliptical distribu- 
tion of soft X-ray emission with radii of 1.3' x 0.75' 
or 19.4 X 11.25 pc. The shell has typical values for 
[S ii]/Hq! of ~ 0.6, as opposed to maximum values 
~ 0.2 elsewhere. The bright Ha emission to the 
east of SNR0449-6921 comes from the Hii region 
N79 and is not directly related to our study. See 
Figures El and H 

N79 corresponds to an area of diffuse soft X- 
ray emission which covers half the XMM field of 
view. We defined a circular source region cover- 
ing the SNR and then divided the field of view 
into two background regions corresponding to the 
area covered by the diffuse X-ray emission and the 
empty area. The spectra from the two background 
regions were fit simultaneously with independent 
normalizations but other model parameters tied 
together. In the region with diffuse X-rays we 
included an additional absorbed thermal compo- 



nent. 

SNR0449-6921 shows a clear bow-shaped ex- 
pansion pattern in the echelle spectrum. It 
is difficult to discern a complete corresponding 
blueshifted pattern. It may be obscured by the 
broad band of emission from the surrounding H ii 
region if SNR0449-6921 is moving away along our 
line of sight with respect to the systemic veloc- 
ity. Some blueshifted emission is seen along the 
right side of the spectrum. Taking this to be the 
front edge of the shell, we measured a velocity of 
±70 km s-^ 

We found a best-fit temperature for the local 
hot bubble of kT = 0.11 keV, and for the Galactic 
halo of kT — 0.27 keV. In the region with the dif- 
fuse soft X-rays we found an additional absorbed 
thermal component with kT = 0.19 keV. Best-fit 
hydrogen column density is 3.7x10^^ cm~^. Full 
details of the background fit are found in Table [H 

Unfortunately all but 7 ks of this observation 
was lost to flaring events. This combined with 
the fact that SNR0449-6921 is small and faint re- 
sulted in too few source counts to constrain the 
SNR properties, even with the information pro- 
vided by the background fits. However, by defin- 
ing emission above 1.5 keV as hard, we can esti- 
mate a hardness ratio of -0.97. Based on the clear 
detection of soft X-ray emission confined within 
a well-defined shell of enhanced [Sil]/Ha optical 
emission, we identify this candidate as a bona fide 
SNR. 

3.2. SNR0506 6542 

SNR0506-6542 (DEM L 72, FHW95 B0505- 
6548) shows a clear shell structure in the opti- 
cal data with [Sii]/Ha between 0.5 and 0.7. It 
is isolated from any other diffuse optical or X- 
ray emission. It is roughly circular along the 
north edge, but flattened along the eastern side. 
The southwestern edge shows a bright fiat region 
with complex filamentary structure. There ap- 
pears to be a gap in the shell along the south- 
ern edge but closer investigation shows faint op- 
tical filaments extending southwestward and re- 
connecting with the fiat region mentioned above. 
Complex filamentary structure is visible inside the 
shell. We also see a filament extending northeast 
from the remnant. The fainter filaments are most 
clearly visible in [O ill] images. The optical shell 
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and X-ray emission cover 7.3' or 110 pc along a 
northwest-southeast hne, and 6.7' or 100 pc along 
a northeast-southwest line. A northeastern spur 
extends a projected distance of 1.2' or 18 pc east 
of the main shell. For the purposes of our calcu- 
lations, we assign a definite size to the remnant 
by fitting an ellipse with radii of 2'.1 and 4.'0, or 
40 pc and 60 pc, to the shell. This ellipse encom- 
passes the main body of the shell, along with the 
northeastern spur and the faint filament bounding 
the southern edge. The shell is filled with diffuse 
soft X-ray emission but it is not as evenly dis- 
tributed as in SNR0449-6921. Several local max- 
ima are observed. The X-ray emission seems to 
divide itself into three regions, one in the north- 
west corner, a second within the southern bulge, 
and a third along the eastern side, fully contained 
by the northeastern spur. The emission from this 
last region seems to be generally fainter than that 
from the other two regions. See Figures [5] and [HI 

We detected emission across the face of SNR0506- 
6542 in the echelle spectra, although the pattern 
is not smooth. Rather it shows several areas of 
convergence toward or deviation from systemic ve- 
locity. The systemic component is distinguished 
from the telluric OH line. Maximum inferred shell 
velocities are ± 90 km s~^. 

In the XMM data, we defined a circular source 
region covering the SNR and extracted the remain- 
der of the FOV as the background. After light 
curve filtering we were left with 9 ks of integra- 
tion time. We found a best fit temperature for 
the local hot bubble as kT = 0.09 keV and for the 
Galactic halo as kT = 0.2 keV, consistent with the 
results from the SNR0449-6921 background. The 
best fit column density is 3.2x10^^ cm^^. Using 
these background parameters, we fit our source re- 
gion with an additional thermal plasma and found 
a gas temperature of kT — 0.17 keV. 

3.3. SNR0537-6628 

SNR0537-6628 (DEM L 256) also lies projected 
near the Hii regions DEM L 251, 253, and 264. 
Optical filaments take the form of two nested 
shells. The inner shell is complete and has a radius 
of 0.8' or 11.6 pc. The outer shell covers roughly 
180° on the southeast half of the remnant with 
a radius of 1.6' or 23 pc along its western extent 
and 1.3' or 20 pc in the perpendicular direction. 
Two clumps of Ha emission appear on the west- 



ern side of the remnant. One clump lies along the 
inner shell, and the second is roughly co-circular 
with the outer shell. Our [S ii]/Hq; ratio maps re- 
veal that these clumps have ratios of < 0.4, similar 
to the nearby Hii regions and in contrast to the 
values > 0.7 found in the shell filaments. Thus 
it seems that these clumps are not actually associ- 
ated with the SNR. The X-ray emission is confined 
to the eastern half of the remnant corresponding 
to the extent of the outer Ha shell and appears to 
be contained primarily between the two filaments. 
See Figures [7] and [HI 

This SNR shows clear blue- and redshifted ex- 
pansion patterns across the face, although the red- 
shifted emission is considerably stronger than the 
blueshift counterpart. A wide band of Ha emission 
is also visible from the surrounding H ii region. 
The expansion pattern has two distinct sections, 
with distinct central velocities but equal expansion 
velocities of ± 55 km s~^. The break between the 
two patterns corresponds to the point where the 
echelle slit crosses the inner optical shell. Strong 
Ha emission is seen converging to the systemic 
velocity at the location of the outer half-shell. A 
corresponding convergence on the other end of the 
slit is much fainter, but this end of the slit does 
not intersect any part of the outer half-shell. 

There is diffuse soft X-ray emission around the 
edge of the XMM field of view in the direction 
of DEM L 256. In this exposure, 20 ks of inte- 
gration time were useable. We extracted a back- 
ground annulus around the SNR that is small 
enough to avoid the emission associated with the 
Hii regions. Here we found best-fit temperatures 
of kT = 0.17 keV for the local hot bubble and 
kT — 0.22 keV for the Galactic halo, again consis- 
tent with the above findings. Column density is fit 
as 5x10^^ cm"^. Fitting the center of the annulus 
with this background plus a source plasma gave a 
gas temperature of kT = 0.22 keV. 

4. Discussion 

The [Sii]/Ha ratio technique has been used 
before to discover SNRs which were too faint 
to be detected i n pre vious surveys, for e x ample , 



Williams et all (|2004D and IStupar et al.l (|2007f ) 



This technique when applied to the detailed and 
deep MCELS images turned up yet more SNR can- 
didates which had been missed in previous optical 
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studiefl Using the thermal source models fitted 
to the X-ray spectra, we can derive a flux which 
can be converted to an estimated ROSAT High 
Resolution Imager count rate using the PIMMS 
tool from HEASARC. We found that ROSAT 
would have collected roughly 1.3 x 10"'* counts 
per second per square arcminute from SNR0506- 
6542 and 3 x lO'^ from SNR0537-6628 over its 
primary energy band of 0.2-2.0 keV. We were not 
able to calculate a flux for SNR0449-6921 since 
its spectrum was of too poor quality to achieve 
a satisfactory fit. However we can set an up- 
per limit on its fiux by assuming that all emis- 
sion from the source region is due to a single ab- 
sorbed thermal plasma. In this way we found a 
flux of 5.0 X 10"^'' erg s~^ cm~^ corresponding to 
2.1 X 10"'^ counts per second per square arcminute 
for ROSAT. We found that the two remnants pre- 
sented here for which we have good flux measure- 
ments are indeed fainter than all ROSAT detected 



LMC remnants from the catalog of lWilliams et al 
()l999l) . Our upper limit for SNR0449-6921 is 



higher than three of the ROSAT SNRs, but among 
those, one benefitted from a lengthy 109 ks ex- 
posure, and the other two detections are of too 
poor quality to allow a discussion in that refer- 
ence. We also note that SNR0537-6628 coincides 
with HP 344, although t hat object was not identi - 
fied as a SNR candidate (|Haberl fc Pietschlll999l ). 
Again using PIMMS to convert our measured flux 
to a ROSAT PSPC count rate over the range 0.1- 
2.4 keV we found an expected rate of 1.6 x 10~^ 
counts per second, which compares favorably with 
the r ate of (1.4 ± 0.2) x 10~ ^ counts per second 
from lHaberl fc PietschI (|l999l ). 

Our technique can not be expected to identify 
the complete p opula t ion of fain t SNRs, however 
As dis cussed in IChul (|l997l) and IChu fc Mac Lowl 
(|l990l ). SNRs in superbubbles cleared out by OB 
associations lack the optical shell signature. Sev- 
eral SNRs have been discovered within LMC su- 
perbubbles (SBs) solely by their X-ray emission 
in ROSAT mosaics. Based on the success of the 
present project in identifying previously unknown 
faint remnants, it seems likely that more such rem- 
nants exist, but may be lacking in optical signa- 
tures. A deep X-ray survey of SBs with modern 



instruments would be necessary to flll out the pop- 
ulation of LMC remnants. Since the spectra of 
SNRs soften as they age, instrumentation with en- 
hanced sensitivity to low energy photons would be 
especially beneficial. 

Superbubbles in the LMC are often in excess of 
100 pc, so the comparatively small sizes of these 
objects, as well as the lack of interior OB associa- 
tions, argues for their formations as single SNRs, 
rather than as SBs. Likewise, SBs tend to have 
fairly slow expansion velocities, close to those of 
the surrounding ISM. The derived expansion ve- 
locities are on the low side for SNRs, but are higher 
than one would expect for a SB toun ne et al.l 
200 1[ ). Only SNR0506-6542 is of comparable size 



to these SBs but has a higher expansion velocity, 
which argues for an SNR origin. SNR0449-6921 
and SNR0537-6628 are on the high side of velocity 
for SBs, are notably smaller than typical SBs, and 
do not have any known OB associations, which 
suggests that they also are likely SNRs. 

We can estimate the density of the ambient ISM 
by summing the masses of the warm shell and hot 
gas and dividing by the volume of the remnant. 
Technically one should also subtract the mass of 
the progenitor star, but in all three cases here the 
SNR masses are now much greater than a typ- 
ical progenitor mass of ^IOMq. Moreover, this 
ISM density estimate will be strongly affected by 
the same approximations mentioned in §2 with re- 
spect to calculating the shell mass, so these are 
at best a rough estimate. For SNR0449-6921, we 
were unable to fit a model to the source plasma 
and so we do not have information on the density 
of the hot gas. In this case, we included only the 
mass of the shell and calculated a lower limit of 
^ 4 cm^"^ on the ambient density. Even this lower 
limit is significantly higher than the densities of 0.6 
and 1.4 cm~'^ which we estimated for SNR0506- 
6542 and SNR0537-6628 respectively. This is con- 
sistent with the observation that SNR0449-6921 
seems to be embedded in an H ii region. 

An additional motivation for studying SNR0449- 
6921 was its close projection to the ring nebula 
around the Wolf-Rayet star Br2. This ring neb- 
ula is unusual in its bright nebular HeiiA4686 
line emission, indica ting very high excitation 



^See for example Williams et al. 2010 (in preparation) and 
references therein, online at ^http:/ /hoth.ccssc.org/mcsnr/] 



([Garnett fc Chul 119941 ) . Two other examples of 
He ll-emitting nebulae are known in the LMC: 
the Hen nebula around the X-ray binary LMC 
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X-1 (jPakull fc Angibaultl Il986l ) and N44C which 
has been s uKgested to be ionized by the transient 



LMC X-5 (IPaknll fc Motchlll989^ . Later studies 
demonstrated that the He ii and Ha hne widths in 
the ring nebula are too narrow to be caused by dy- 
namical interacti on with the expanding SNR shell 
(|Chu et alJll999[ ). Jowever, the possibility of ex- 



citation by an X -ray source within the ring (e.g., 
Chu et allbood ) or within the SNR had not been 
ruled out. Our obseryations found no significant 
compact X-ray emission at the location or Br2. 
Furthermore, we see a smooth distribution of soft 
emission within the SNR shell with no indications 
of any compact or hard source. 

The very irregular shell morphology of SNR0506- 
6542 makes it difficult to assign it an exact size. 
Measured at the extreme extents of the optical 
shell, it has radii of 50 x 55 pc which would 
make SNR0506-6542 the largest known SNR 
i n the LMC, l arger even than SNR 0450-70.9 
(jWilliams et al.l 12004) . However, for the purposes 
of our quantitative estimates, we take a shell size 
of 40 X 60 pc, as mentioned above. The bright- 
est X-ray emission is seen in the western half of 
the SNR. Here the optical shell flattens and be- 
comes thicker and more filamentary. Many optical 
clumps are also seen in this half of the shell. This 
is consistent with the SNR encountering a region 
of higher ISM density which is being shocked by 
the blast wave. As this material passes through 
the shock, it evaporates producing diffuse X-ray 
emission. The faint optical filament along the ex- 
treme eastern edge of the remnant coincides with 
a line of X-ray emission. This shock front may 
still be capable of heating the gas it encounters to 
X-ray emitting temperatures. Elsewhere along the 
northern rim we see less emission along the shell 
and, in some areas, a gap between the optical shell 
and interior X-ray emission, suggesting that the 
shock front here is no longer powerful enough to 
fully ionize new gas and that previously heated 
gas has now cooled to below detectable levels. 

The expansion velocity of 90 km s^^ implies an 
age of ~150 kyr. Solving for E/n as described in 
§2 gives E/n — 2.6 x 10^^ erg cm^, in excellent 
agreement with the values calculated from the Ha 
surface brightness and X-ray spectrum as listed in 
Table [2] However, such slow- moving shocks would 
not be capable of heating new gas to X-ray emit- 
ting temperatures. Post-shock velocities of '--^140 



km s~^ are needed to heat gas to a sufficient tem- 
perature to radiate strongly at 0.2 keV. The evi- 
dence for newly shock-heated gas seen in the X- 
ray images may indicate that for at least some 
portions of the SNR front, the velocity is well in 
excess of 90 km s~^. If we instead consider the 
expansion velocity to be 140 km s~^, we find an 
age of 95 kyr. This is an issue seen in other rem- 
nants, e.g. N86 and NllL. It may be that the op- 
tical material with the greatest Doppler shifts has 
too low surface brightness for detection in echelle 
spectra. In any case, the observed gas velocities 
are well below those expected at the end of the 
Sedov- Taylor phase, and thus the age of 150 kyr 
is likely an overestimate. For comparison, an SNR 
of this radius but with gas velocity corresponding 
to the end of the Sedov- Taylor phase would give 
an age of '^80 kyr. 

When applying the simulated spectrum, we 
found a best fit at an age of 80 kyr. However, the 
simulations seem to underprcdict the prominent 
emission lines around 0.6 keV. This could indi- 
cate that we have underestimated the abundances 
throughout our fitting process. Returning to the 
background model and allowing the local hot bub- 
ble abundance to float, we found a best fit when 
ZhHB = 1-9. Applying the background model to 
the source spectrum we again found a best fit to 
the SNR emission when the source plasma abun- 
dances are set to 0.5. Plugging in the simulated 
spectrum we found that an age of 100 kyr is now 
preferred, in better agreement with the analyti- 
cally derived values. However, one might expect 
that such an old remnant would be heavily dom- 
inated by the swept-up medium, and therefore 
show abundances close to those of the LMC. Fur- 
thermore, such an old remnant will have evolved 
past the Sedov- Taylor phase which ends at veloc- 
ities 

It is worth pointing out that the physical 
parameters which we derived for this SNR are 
broadly similar to those found by Williams et al.l 
()2004l) for SNR 0450-70.9. Prior to that study. 



no interior X-ray emission had ever been detected 
from a remnant of this size. Moreover, no remnant 
of this size is known within our galaxy. Thus, a 
detailed study of two such advanced SNRs repre- 
sents a significant step forward in our understand- 
ing of the long-term fate of these objects. We also 
note that Phot is significantly greater than PghoU 
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for this remnant, suggesting some pressure-driven 
expansion even at these late ages. 

X-ray emission from SNR0537-6628 appears to 
trace the Ha filaments, again suggesting that this 
gas is recently shock heated. The appearance of 
two distinct shells along with two expansion pat- 
terns in the echelle spectrum which correlate with 
the optical shells suggests a bilobed structure for 
the SNR. The near edge of the SNR seems to be 
encountering denser material as evidenced by the 
X-ray emission corresponding to that shock wave 
and the stronger Ha emission converging to sys- 
temic velocity on that edge. 

The observed expansion velocity of 55 km 
is quite low and gives an age of 125 kyr. The 
corresponding ratio E/n = 6 x 10^^ erg cm^ 
docs not agree as well with the observed value of 
1.6 X 10^° erg cm"^. However, the thermal energy 
of the hot gas was estimated assuming filling of an 
entire ellipsoidal shell, whereas we see X-ray emis- 
sion from only the eastern half. Furthermore, if 
this is a bilobed structure which is not aligned with 
the line of sight, we are only measuring a projected 
component of the true expansion velocity. Both of 
these considerations would bring the two numbers 
into closer agreement. We calculate that the view- 
ing angle of the bilobed expansion would need to 
be 40° in order for the two numbers to agree. 
We also note that if the observed elliptical region 
of X-ray emission were actually a circular disc, a 
viewing angle of ~ 55° would cause it to appear as 
it does. This viewing angle would imply a true ex- 
pansion velocity of 97 km s~^ which corresponds 
to an age of 70 kyr. However, if the remnant is 
indeed cncoimtcring denser material, it may have 
undergone a more abrupt (non-adiabatic) deceler- 
ation, and again, this expansion velocity is outside 
the expected range of the Sedov- Taylor phase, net- 
ting a younger age for a given observed radius and 
expansion velocity. 

The imcertainty in the properties of this SNR 
makes it difficult to choose accurate parameters 
for the simulations. Moreover, our simulations 
arc spherically symmetric, which is clearly not 
the case for this remnant. Simulations run over 
a range of reasonable parameters all result in fits 
with ages < 50 kyr, lending support to the above 
consideration. 



5. Conclusions 

We have presented a multi-wavelength study of 
three recently discovered SNRs in the LMC, fea- 
turing new XMM-Newton observations. We fit 
background spectra from each observation with a 
detailed model including soft proton contamina- 
tion, instrumental lines, the local hot bubble, and 
extragalactic components. These fits gave consis- 
tent values for the background parameters across 
the three observations. We then fit source spectra 
with thermal plasma models, and compare the fits 
to simulated spectra. We analytically estimated 
the physical properties of the SNRs. 

SNR0449-6921 is confirmed as an SNR due to 
its [SIT] /Ha ratio, velocity expansion pattern, and 
soft X-ray emission, although we cannot fully char- 
acterise that emission. SNR0506-6542 appears 
from our findings to be a large, relatively old SNR, 
similar to SNR 0450—709. It has fairly disorga- 
nized expansion and X-ray emission that appears 
to be due to a combination of recently shocked 
material and older "fossil" radiation. SNR0537- 
6628 appears to be encountering higher density 
material to one side and shows signs of a bipolar 
expansion geometry. Although we were able to ob- 
tain good fits to its spectrum, the age is difficult 
to determine precisely due to uncertainties in its 
geometry and environment. 
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data and folded model 




0.5 1 2 

channel energy (keV) 



Fig. 1. — The background spectrum for SNR0506- 
6542 with individual components plotted and la- 
beled. 
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channel energy (keV 



Fig. 2. — Background-subtracted source spectrum 
for SNR0506-6542 fitted with our source model. 
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4:50:00.0 



40.0 



20.0 



Right ascension 

Fig. 3.— Three-color image of SNR0449-6921 
with X-ray contours overlaid. Red is Ha, green 
is [Sii], and blue is [Oiii]. Contours are at lev- 
els of 2 counts per second per square degree in 
the energy band 0.35-0.85 keV. The dashed ellipse 
corresponds to the SNR extent considered in the 
calculations of the X-ray and optical shell proper- 
ties. The line corresponds to the radial slice shown 
in Figure 01 
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Fig. 4.— Radial profile of SNR0449-6921. The 
solid line represents Ha and the dashed line rep- 
resents 0.35-0.85 keV smoothed X-rays. X-axis 
units are parsccs. The origin corresponds to 04 
49 22.74, -69 20 27.28. The slice was chosen so as 
to avoid the diffuse emission from the nearby H ii 
region. 



c 
o 
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40.0 20.0 5:06:00.0 40.0 05:20.0 

Right ascension 



Fig. 5.— Three-color image of SNR0506-6542 
with X-ray contours overlaid. Red is Ha, green 
is [Sii], and blue is [Oiii]. Contours are at lev- 
els of 2 counts per second per square degree in 
the energy band 0.35-0.85 keV. The dashed ellipse 
corresponds to the SNR extent considered in the 
calculations of the X-ray and optical shell proper- 
ties. The line corresponds to the radial slice shown 
in Figure ini 
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Fig. 6.— Radial profile of SNR0506-6542. The 
solid line represents Ha and the dashed line rep- 
resents 0.35-0.85 keV smoothed X-rays. X-axis 
units arc parsecs. The origin corresponds to 05 
06 01.87, -65 41 29.43. The northeastern spur is 
visible near -45. 



10.0 38:00.0 50.0 40.0 30.0 20.0 10.0 5:37:00.0 50.0 36 

Fig. 7.— Three-color image of SNR0537-6628 
with X-ray contours overlaid. Red is Ha, green 
is [S ii] , and blue is [O ill] . Contours are at lev- 
els of 2 counts per second per square degree in 
the energy band 0.35-0.85 keV. The dashed ellipse 
corresponds to the SNR extent considered in the 
calculations of the X-ray and optical shell proper- 
ties. The line corresponds to the radial slice shown 
in Figure [SI 
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Fig. 8.— Radial profile of SNR0537-6628. The 

solid line represents Ha and the dashed line rep- 
resents 0.35-0.85 keV smoothed X-rays. X-axis 
units are parsecs. The origin corresponds to 05 37 
28.49, -66 27 41.47. The shoe was chosen to lie 
along the western extent of the outer shell. 



Table 1 
XMM-Newton Observation Log 



ObsID 


RA 


Dec. 


Start Time 


Exp. Time (ks) 


205260101 


04h49ni40.00s 


-69°21'49.0" 


2004-02-24 16:19:07 


18.7 


205260201 


05h06m05.00s 


-65°41'30.0" 


2004-01-16 18:06:28 


16.9 


301410601 


05h37m32.00s 


-66°28'26.6" 


2005-10-02 09:33:11 


21.8 
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Table 2 

Physical Properties of SNRs. The X-ray fit parameters are from an APEC model. 





SNR0449 6921 


SNR0506-6542 


SNR0537 -6628 


R (pc) 


19x11 


60x40 


24x20 


SB (10-15) 


0.9 


0.5 


0.4 


Wexp (km s-i) 


70 


90 


55 


"e,shell (cm-3) 


5.5(2) 


2.99(8) 


3.30(4) 




1.4 X 10^9 


1.7 X 10^0 


3.2 X 10^9 




820(30) 


5.6(1) x 10-'' 


1.14(l)xl03 


-K' (erg) 


4(1) X 10^9 


4(1) X 10^° 


3(1) X 10*9 


-Psheii (dyne cm- 2) 


1.51(5) X 10-11 


8.2(2) X 10-12 


9.1(1) X 10-12 


/cT (kcV) 




0.17(1) 


0.22(5) 


Normalization (IQ-^'* cm-^) 




0.004(2) 


0.0009(9) 


V (cm3) 




1.3(3) X lO'^i 


1.0(1) X 10^0 


J^c.hot (cm-3) 




0.11(4) 


0.2(2) 
2(2)xl02 


Mhot (Mo) 




1.3(9) X 10^ 


-E'th (erg) 




1.1(8) X 10^1 


2(2) X 10^° 


-Phot (dyne cm-^) 




5(2) X 10-11 


1(1) X 10-1° 


"ISM (cm-3) 


> 4.0 


0.5 1(1) 


1.40(2) 
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Table 3 

Background Spectral Parameters 





SNR0449-6921 
Hll region 


SNR0449-6921 
clear region 


SNR0506-6542 


SNR0537-6628 


Soft Proton Index 


1.9(6) 


1.9(6) 


1.0(1) 


1.39(7) 


Soft Proton Norm, 
(cts keV~^ cm~^ s~-^) 


0.02(2) 


0.003(3) 


0.04(1) 


0.048(5) 


Local Bubble Temp. (keV) 


0.11(1) 


0.11(1) 


0.09(1) 


0.17(4) 


Local Bubble Norm. (IQ-^^ cm-^) 


4(4) X 10-4 


7(7) X 10-4 


0.002(2) 


3(1) X 10-4 


A/-H (10^2 cm-2) 


0.37(7) 


0.37(7) 


0.32(7) 


0.5(1) 


Extragal. Thermal Temp. (keV) 

Extragal. Thermal Norm. 
(10-" cm-5) 


0.27(4) 
0.19(4) 

7(7) X 10-4 
0.02(1) 


0.27(4) 
0.001(1) 


0.20(2) 
0.003(2) 


0.22(2) 
0.005(3) 


Extragal. Norm. 

(10-^ cts keV-i cm-2 s"!) 


2(2) 


2.1(7) 


3(2) 


0.9(9) 
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